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The influence of helicity on the stability of scaling regimes, on the effective diffusivity, and on 
the anomalous scaling of structure functions of a passive scalar advected by a Gaussian solenoidal 
velocity field with finite correlation time is investigated by the field theoretic renormalization group 
and operator product expansion within the two-loop approximation. The influence of helicity on 
the scaling regimes is discussed and shown in the plane of exponents e — rj, where e characterizes 
the energy spectrum of the velocity fleld in the inertial range E oc k^~^^, and rj is related to 
the correlation time at the wave number k, which is scaled as k~^^^ . The restrictions given by 
nonzero helicity on the regions with stable fixed points that correspond to the scaling regimes are 
analyzed in detail. The dependence of the effective diffusivity on the helicity parameter is discussed. 
The anomalous exponents of the structure functions of the passive scalar field which define their 
anomalous scaling are calculated and it is shown that, although the separate composite operators 
which define them strongly depend on the helicity parameter, the resulting two-loop contributions to 
the critical dimensions of the structure functions are independent of helicity. Details of calculations 
are shown. 
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I. INTRODUCTION 

During the last decade much attention has been paid 
to the inertial range of fully developed turbulence, which 
contains wave numbers larger than those that pump 
the energy into the system and smaller than those that 
are related to the dissipation processes [ll, [1]. Ground- 
ing for the inertial range turbulence was created in the 
well known Kolmogorov-Obukhov (KO) phenomenologi- 
cal theory (see, e.g., [H, H, One of the main problems 
in the modern theory of fully developed turbulence is 
to verify the validity of the basic principles of KO the- 
ory and their consequences within the framework of a 
microscopic model. Recent experimental and theoretical 
studies indicate possible deviations from the celebrated 
Kolmogorov scaling exponents. The scaling behavior of 
velocity fluctuations with exponents whose values are dif- 
ferent from Kolmogorov ones is called anomalous and 
usually is associated with the intermittency phenomenon. 
In turbulence this phenomenon is believed to be related 
to strong fluctuations of the energy flux which there- 
fore leads to the deviations from the predictions of the 
aforementioned KO theory. Such deviations, referred to 
as anomalous or nondimensional scaling, manifest them- 
selves in a singular dependence of the correlation or struc- 
ture functions on the distances and the integral (exter- 
nal) turbulence scale L. The corresponding exponents 



are certain nontrivial and nonlinear functions of the or- 
der of the correlation function, the phenomenon referred 
to as "multiscaling." Even though great progress in the 
understanding of intermittency and anomalous scaling in 
turbulence was achieved as a result of intensive studies, 
their investigation in fully developed turbulence still re- 
mains a major theoretical problem. 

Although the theoretical description of fluid turbulence 
on the basis of first principles, i.e., on the stochastic 
Navier-Stokes equation [1], remains essentially an open 
problem, considerable progress has been achieved in un- 
derstanding simplified model systems that share some 
important properties with the real problem: shell mod- 
els the stochastic Burgers equation and passive 
advection by random "synthetic" velocity fields [3]. 

A crucial role in these studies is played by mod- 
els of the advected passive scalar field A simple 
model of a passive scalar quantity advected by a random 
Gaussian velocity field, white in time and self-similar 
in space (the latter property mimics some features of a 
real turbulent velocity ensemble), the so-called Kraich- 
nan rapid-change model [9], is an example. The interest 
in these models is based on two important facts: first, as 
shown by both natural and numerical experimental in- 
vestigations, the deviations from the predictions of the 
classical Kolmogorov-Obukhov phenomenological theory 
[E H, 0, are even more strongly displayed for a pas- 
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sively advected scalar field than for the velocity field it- 
self (see, e.g., [nl, [H [H, H [S [3 and references cited 
therein) , and second, the problem of passive advection is 
much easier to consider from a theoretical point of view. 
In these studies the anomalous scaling was established on 
the basis of a microscopic model \l7^ , and corresponding 
anomalous exponents were calculated within controlled 
approximations [H, [H, 2Q| (see also the review Q and 
references therein). 

The greatest stimulation to study the simple models 
of passive advection not only of scalar fields but also of 
vector fields (e.g., a weak magnetic field) is related to 
the fact that even simplified models with given Gaussian 
statistics of a so-called synthetic velocity field describes 
a lot of features of the anomalous behavior of genuine 
turbulent transport of some quantities (such as heat or 
mass) observed in experiments (see , e.g., Refs. fol. ITsl. IT^. 
m, III [13, [II, III S [111, 111 til and references cited 
therein) . 

An effective method for investigation of a self-similar 
scaling behavior is the renormalization group (RG) tech- 
nique ^25,, 2^ ■ It was widely used in the theory of crit- 
ical phenomena to explain the origin of the critical scaling 
and also to calculate corresponding universal quantities 
(e.g., critical dimensions). This method can also be di- 
rectly used in the theory of turbulence [H [13, HI S HI , 
as well as in related models like the simpler stochas- 
tic problem of a passive scalar advected by a prescribed 
stochastic flow. In what follows we use the conventional 
("quantum field theory" or field theoretic) RG which is 
based on the standard renormalization procedure, i.e., on 
elimination of the ultraviolet (uv) divergences. 

In Refs. [3l| the field theoretic RG and operator- 
product expansion (OPE) were used in the systematic in- 
vestigation of the rapid-change model. It was shown that 
within the field theoretic approach the anomalous scal- 
ing is related to the very existence of so-called dangerous 
composite operators with negative critical dimensions in 
the OPE (see, e.g., [25, 30] for details). In subsequent 
papers [32| the anomalous exponents of the model were 
calculated within the e expansion to order (three-loop 
approximation). Here e is a parameter that describes a 
given equal-time pair correlation function of the veloc- 
ity field (see later section). Important advantages of the 
RG approach are its universality and calculational effi- 
ciency: a regular systematic perturbation expansion for 
the anomalous exponents was constructed, similar to the 
well-known e expansion in the theory of phase transitions. 

Afterward, various generalized descendants of the 
Kraichnan model, namely, models with inclusion of large 
and small scale anisotropy 33] , compressibility 13411 , and 
finite correlation time of the velocity field [35l 1361 ] were 
studied by the field theoretic approach. Moreover, ad- 
vection of a passive vector field by the Gaussian self- 
similar velocity field (with and without large and small 
scale anisotropy, pressure, compressibility, and finite cor- 
relation time) has also been investigated and all possible 
asymptotic scaling regimes and crossovers among them 



have been classified [37*]. The general conclusion is that 
the anomalous scaling, which is the most important fea- 
ture of the Kraichnan rapid-change model, remains valid 
for all generalized models. 

Let us describe briefiy the solution of the problem in 
the framework of the field theoretic approach (see, e.g., 
Refs. fl^, [1^, [3^] for more details). It can be divided 
into two main stages. In the first stage the multiplica- 
tive renormalizability of the corresponding field theoretic 
model is demonstrated and the differential RG equations 
for its correlation functions are obtained. The asymptotic 
behavior of the latter on their ultraviolet argument {r/£) 
for r ^ £ and any fixed {r/L) is given by infrared stable 
fixed points of those equations. Here i and L are the 
inner (ultraviolet) and the outer (infrared) scales. The 
behavior involves some "scaling functions" of the infrared 
argument (r/L), whose form is not determined by the RG 
equations. In the second stage, their behavior at r <C L 
is found from the OPE within the framework of the gen- 
eral solution of the RG equations. There, the crucial role 
is played by the critical dimensions of various composite 
operators, which give rise to an infinite family of indepen- 
dent scaling exponents as mentioned above (and hence to 
multiscaling). Of course, both these stages (and thus the 
phenomenon of multiscaling) have long been known in 
the RG theory of critical behavior. The distinguishing 
feature specific to models of turbulence is the existence 
of composite operators with the aforementioned negative 
critical dimensions. Their contributions to the OPE di- 
verge at {r/L) — > 0. In models of critical phenomena, 
nontrivial composite operators always have strictly posi- 
tive dimensions, so that they only determine corrections 
[vanishing for {r/L) —> 0] to the leading terms [finite for 
{r/L) 0] in the scaling functions. 

In Ref. [11 the problem of a passive scalar advected 
by a Gaussian self-similar velocity field with finite cor- 
relation time [H was studied by the field theoretic RG 
method. There, a systematic study of the possible scaling 
regimes and anomalous behavior was presented at one- 
loop level. The two-loop corrections to the anomalous 
exponents were obtained in Ref. 39]. It was shown that 
the anomalous exponents are nonuniversal as a result of 
their dependence on a dimensionless parameter, the ratio 
of the velocity correlation time and the turnover time of 
the scalar field. 

In what follows, we shall continue with the investiga- 
tion of this model from the point of view of the influence 
of helicity (spatial parity violation) on the scaling regimes 
and anomalous exponents within the two-loop approxi- 
mation. 

Helicity is defined as the scalar product of velocity and 
vorticity and its nonzero value expresses mirror symme- 
try breaking of the turbulent fiow. It plays a significant 
role in the processes of magnetic field generation in an 
electrically conductive fiuid [il-lil and represents one 
of the most important characteristics of large-scale mo- 
tions as well [47] -[11. The presence of helicity is observed 
in various natural (like large air vortices in the atmo- 
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sphere) and technical flows [i^ [U HI] . Despite this fact 
the role of the helicity in hydrodynamical turbulence is 
not completely clarified up to now. 

The Navier-Stokes equations conserve kinetic energy 
and helicity in the inviscid limit. The presence of two 
quadratic invariants leads to the possibility of appear- 
ance of double cascade. This means that cascades of 
energy and helicity take place in different ranges of wave 
numbers analogously to the two-dimensional turbulence 
and/or the helicity cascade appearing concurrently to the 
energy cascade in the direction of small scales [H, HSl- 
In particular, the helicity cascade is closely connected 
with the existence of the exact relation between triple 
and double correlations of velocity known as the 2/15 law 
analogously to the 4/5 Kolmogorov law [5^ . Correspond- 
ing to [1^ the aforementioned scenarios of turbulent cas- 
cades differ from each other by spectral scaling. The- 
oretical arguments given by Kraichnan [56] and results 
of numerical calculations of the Navier-Stokes equations 
(stI . [ssI . [sol support the scenario of concurrent cascades. 
The appearance of helicity in turbulent systems leads to 
the constraint of the nonlinear cascade to small scales. 
This phenomenon was first demonstrated by Kraichnan 
[56j within the modeling problem of statistically equilib- 
rium spectra and later in numerical experiments. 

Turbulent viscosity and diffusivity, which characterize 
the influence of small-scale motions on heat and momen- 
tum transport, are basic quantities investigated in the 
theoretic and applied models. The constraint of the di- 
rect energy cascade in helical turbulence has to be ac- 
companied by a decrease of turbulent viscosity. How- 
ever, no influence of helicity on turbulent viscosity was 
found in some works [60L [6ll ] . A similar situation is ob- 
served for the turbulent diffusivity in helical turbulence. 
Although the modeling calculations demonstrate intensi- 
fication of turbulent transfer in the presence of helicity 
[6^ . [63t direct calculation of the diffusivity does not con- 
firm this effect [H, [H, HI). Helicity is a pseudoscalar 
quantity, hence, it can be easily understood that its in- 
fluence appears only in quadratic and higher terms of 
perturbation theory or in combination with other pseu- 
doscalar quantities (e.g., large-scale helicity). Really, si- 
multaneous consideration of memory effects and second 
order approximation indicates the effective influence of 
helicity on turbulent viscosity [H, 113] and turbulent dif- 
fusivity [H, [H, [zS] already in the limits of small and 
infinite correlation time. 

Helicity, as we shall see below, does not affect known 
results in the one-loop approximation and, therefore, it 
is necessary to turn to the second order (two-loop) ap- 
proximation to be able to analyze possible consequences. 
It is also important to say that in the framework of the 
classical Kraichnan model, i.e., a model of passive ad- 
vection by a Gaussian velocity field with (5-like correla- 
tions in time, it is not possible to study the infiuence of 
the helicity because all potentially helical diagrams are 
identically equal to zero at all orders in the perturbation 
theory. In this sense, the investigation of the helicity 



in the present model can be consider as the first step 
toward analyzing the helicity in genuine turbulence. In 
fact, it is interesting and important to study the helic- 
ity effects because many turbulence phenomena are di- 
rectly infiuenced by them (like large air vortices in the 
atmosphere). For example, in stochastic magnetic hy- 
drodynamics, which studies the turbulence in electrically 
conducting fiuids, it leads to the nontrivial fact of the 
existence of a so-called turbulent dynamo - the genera- 
tion of a large-scale magnetic field by the energy of the 
turbulent motion [13, M, US, S S, Si, [11] . This is an 
important effect in astrophysics. 

The main result of the paper will be the conclusion 
that helicity does not change the anomalous exponents 
of the single-time structure functions within the two-loop 
approximation although the separate composite opera- 
tors which define them strongly depend on the helicity 
parameter. This result leads to the following interesting 
but nontrivial question: is this result related only to the 
two-loop calculations or does it hold for all orders of per- 
turbation theory? Of course, the answer to this question 
definitely lies at least within the three-loop approxima- 
tion. On the other hand, as will be shown, the effective 
diffusivity rather strongly depends on the helicity param- 
eter. 

The paper is organized as follows. In Sec.|IT]we present 
the definition of the model and introduce the helicity 
to the transverse projector of a given pair correlation 
function of the velocity field. In Sec. lIIII we give the 
field theoretic formulation of the original stochastic prob- 
lem and discuss the corresponding diagrammatic tech- 
nique. In Sec. lIVI we analyze the ultraviolet divergences 
of the model, establish its multiplicative renormalizabil- 
ity, and calculate the renormalization constants in the 
two-loop approximation. In Sec.|V] we analyze possible 
scaling regimes of the model, associated with nontrivial 
and physically acceptable fixed points of the correspond- 
ing RG equations. There are five such regimes, any one 
of which can be realized in dependence on the values of 
the parameters of the model. We discuss the physical 
meaning of these regimes (e.g., some of them correspond 
to zero, finite, or infinite correlation time of the advect- 
ing field) and their regions of stability in the space of the 
model parameters. In Sec. lVIl the two-loop corrections 
to the effective diffusivity are calculated. In Sec. lVIII the 
renormalization of needed composite operators is done 
and their explicit dependence on the helicity parameter 
is shown. In Sec. IVIIll discussion of the results is present. 



II. THE MODEL 

In what follows, we shall consider the advection of a 
passive scalar field 9 = d{x) = 9{t,x.) which is described 
by the following stochastic equation: 

dte + v,d,e = j,oM + f, (1) 
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where dt = d/dt, di = d/dxi, vq is the coefficient of 
molecular difFusivity (hereafter all parameters with a sub- 
script denote bare parameters of the unrenormalized 
theory; see below), A = is the Laplace operator, 
Vi = Vi{x) is the ith component of the divergence- free 
(owing to the incompressibility) velocity field v(a;), and 
j9 = f^{x) is a Gaussian random noise with zero mean 
and correlation function 



{f{x)f{x'))^5{t~t')C{vlL), 



(2) 



where the angular brackets (...) hereafter denote the av- 
erage over the corresponding statistical ensemble. The 
noise maintains the steady state of the system but the 
concrete form of the correlator is not essential. The only 
condition that must be satisfied by the function C{y/L) 
is that it must decrease rapidly for r = |r| 3> L, where 
L denotes an integral scale related to the stirring. In the 
case when C depends not only on the modulus of the vec- 
tor r but also on its direction, it plays the role of a source 
of large-scale anisotropy, whereupon the noise can be re- 
placed by a constant gradient of scalar field. Equation 
dl]) then reads (see, e.g., Ref. (35| ) 



dtO + v^^^9 = voAe - h • V. 



(3) 



Here, 9{x) is the fluctuation part of the total scalar field 
Q{x) = 9{x) -I- h • X, and h is a constant vector that 
determines the distinguished direction. The direct for- 
mulation with a scalar gradient is even more realistic; 
see, e.g. Refs. 0, IH 0, il, SH . 

In real problems the velocity field v(a:) satisfies the 
stochastic Navier- Stokes equation. In spite of this fact, 
in what follows, we shall suppose that the velocity field 
is driven by the simple linear stochastic equation p^ . [ssj 



(4) 



where R = R{x) is a linear operation to be specified 
below and = fi{x) is an external random stirring 
force with zero mean and the correlator 



{n{x)f]{x')) ^ d{^{x-x') 
dujd'^k 



xexp[-i(t-<')+ik(x-x')], (5) 

where k = |k| is the wave number, lo is the frequency, d 
is the dimensionality of the x space (of course, when one 
investigates a system with helicity the dimension of the 
X space must be strictly equal to 3; nevertheless, in what 
follows, we shall retain the d-dimcnsionality of all results 
that are not related to helicity so that we can also study 
the d dependence of the nonhelical case of the model). 
The transition to a helical fiuid corresponds to the giving 
up of conservation of spatial parity, and technically this 
is expressed by the fact that the correlation function is 
specified in the form of a mixture of a true tensor and a 



pseudotensor. In our approach, it is represented by two 
parts of the transverse projector 



pp. = 

2J 



F,,(k)+i?,,(k), 



(6) 



which consists of the nonhelical standard transverse pro- 



jector Py(k) 



kikj/k"^ and Hij(k) — ipsijiki/k, 



which represents the presence of helicity in the fiow. 
Here, Siji is Levi-Civita's completely antisymmetric ten- 
sor of rank 3 [it is equal to 1 or —1 according to whether 
(i, j, I) is an even or odd permutation of (1, 2, 3) and zero 
otherwise] , and the real parameter of helicity, p, charac- 
terizes the amount of helicity. Due to the requirement of 
positive definiteness of the correlation function the abso- 
lute value of p must be in the interval \p\ S (0, 1) (i^list. 
Physically, the nonzero helical part (proportional to p) 
expresses the existence of nonzero correlations (v • rot v) . 

We choose the correlator in Eq. ([5]) to be a 5 func- 
tion in time, which is equivalent to the condition that i?-^ 
is independent of frequency [13|] (see also Refs.fssl. \3^). 
Following [3^ |3^ , we shall work with 



and 



R{k)^uoVo{e +m^Y~'^'^, (8) 



the wave-number representation of R(x). Here, the pos- 
itive amplitude factors 50 and wq play the roles of the 
coupling constants of the model, the analogs of the cou- 
pling constant Aq in the Ao(/?* model of critical behavior 
|24l[25|. In addition, 170 is a formal small parameter of the 
ordinary perturbation theory. The positive exponents e 
and rj [e = 0{rf)\ are small RG expansion parameters, the 
analogs of the parameter e = 4 — d in the AoiyS^ theory. 
Thus, we have a kind of double expansion model in the 
£ — 7] plane around the origin e = 77 = 0. An integral 
scale L = 1/m is introduced to provide infrared (ir) reg- 
ularization. In the limit k ^ m the functions ^ and ^ 
take on simple powerlike forms 



Df{Lj,k)^goiy^k 



R{k) = uoi^ofc^"", (9) 



which will be used in calculations in what follows. The 
needed ir regularization will be given by restrictions on 
the region of integration. 

From Eqs. (|4]), ([5]), and ([9]) we receive the statistics of 
the velocity field v. It obeys a Gaussian distribution with 
zero mean and correlator 

{v,{x)v,{x')) = Dl{x-x') 

(S^^-W^^^'^'^) 
X exp[-iw(i - t') + ik(x - x')],(10) 



with 



3f^4-d-2s-ri 



(11) 
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The correlator (fTTj) is directly related to the energy spec- 
trum via the frequency integral [U, [35| 

/2 
dwb%uj,k)~^k^-^'. (12) 
Wo 

Therefore, the coupling constant go and the exponent 
e describe the equal-time velocity correlator or, equiv- 
alently, the energy spectrum. On the other hand, the 
constant Uq and the second exponent rj are related to the 
frequency w ~ uoVok^~^ [or to the function i?(fc), the 
reciprocal of the correlation time at the wave number 
k] which characterizes the mode k [2l|, [H, [H, [tiI, [73 |. 
Thus, in our notation, the value e = 4/3 corresponds to 
the well-known Kolmogorov five-thirds law for the spa- 
tial statistics of velocity field, and rj — 4/3 corresponds 
to the Kolmogorov frequency. Simple dimensional anal- 
ysis shows that the parameters (charges) go and uq are 
related to the characteristic ultraviolet momentum scale 
A (of the order of the inverse Kolmogorov length) by 

go^A^'+\ uo~A'^. (13) 

In Ref. 13] it was shown that the linear model ^ 
[and therefore also the Gaussian model (|10p . (|lip ] is not 
Galilean invariant and, as a consequence, it does not take 
into account the self-advection of turbulent eddies. As 
a result of these so-called sweeping effects the different 
time correlations of the Eulerian velocity are not self- 
similar and depend strongly on the integral scale; see, 
e.g., Ref. [t^. But, on the other hand, the results pre- 
sented in Ref. [l^ show that the Gaussian model gives a 
reasonable description of the passive advection in the ap- 
propriate frame, where the mean velocity field vanishes. 
One more argument to justify the model pil)l . (fTTj) is 
that, in what follows, we shall be interested in the equal- 
time, Galilean-invariant quantities (structure functions), 
which are not affected by the swee pin g, and therefore, as 
we expect (see, e.g., Refs. [s^ IsgI. l39l|). their absence in 
the Gaussian model pO]) . pT|) is not essential. 

At the end of this section, let us briefly discuss two 
important limits of the considered model pU|) , pT|) . The 
first is the so-called rapid-change model limit when uq — > 
oo and gQ = go/uQ — const, 

&{w,k)^g'„,,,k~''-^'+\ (14) 

and the second is the so-called quenched (time- 
independent or frozen) velocity field limit, which is de- 
fined by Mo — > and gQ = go/ua = const, 

D^ito,k) ^ g'y,7rSiLu)k-^+^-'-, (15) 

which is similar to the well-known models of the random 
walks in random environment with long-range correla- 
tions; see, e.g., Refs. [tI. [75|. 



III. FIELD THEORETIC FORMULATION OF 
THE MODEL 

For completeness of our text in this and the next sec- 
tion we shall present and discuss the principal moments 
of the RG theory of the model defined by Eqs. ©, (fTU)) . 
and (fTTj). 

We start with the reformulation of the stochastic prob- 
lem ©-([S]), according to the well-known general theorem 
(see, e.g., Refs.[2^[25j). into the equivalent field theoretic 
model of the doubled set of fields $ = {9,9' with 
the following action functional: 

5($) = ^Jdh d'^xi dt2 d'^x2 

v'^ih, ^i)d{^ (t 1 , xi ; i2 , X2)w; (i2 , X2 ) (16) 
+ J dt d'^x 9' [-dt9 - V.A9 + v^M - h • v] 

+ J dtd'^x v\ {-dt ~R\v,, 

where is defined in Eq. ^ , 9' and v' are auxiliary 
scalar and vector fields, and summations are implied over 
the vector indices. 

It is standard that the formulation through the action 
functional (fT6|) replaces the statistical averages of random 
quantities in the stochastic problem (I3|)-(l5]) with equiva- 
lent functional averages with weight exp5($). The gen- 
erating functionals of the total Green's functions G(A) 
and connected Green's functions W(A) are then defined 
by the functional integral 

G{A) = e^f-^) ^ e^(*)+^*, (17) 

where A{x) — {A^ , A^ , A"^, A"^ } represents a set of arbi- 
trary sources for the set of fields $, 1?$ ee V9V9'VwW 
denotes the measure of functional integration, and the 
linear form ^4$ is defined as 

A$ = j dx[A\x)9{x)+A^\x)9\x) 

+ A'i{x)v.{x) + Ai {x)v[{x)]. (18) 

Following the arguments in [35| . we can put A^ = 
in Eq. (jlSp and then perform an explicit Gaussian inte- 
gration over the auxiliary vector field v' in Eq. ([17]) as 
a consequence of the fact that, in what follows, we shall 
not be interested in the Green's functions involving the 
field v'. After this integration one is left with the field 
theoretic model described by the functional action 

V, {h , Xi ) [Dj;- {h , Xi ; t2 , X2)]"'wj (t2, X2) (19) 

+ [ dt d'^x 9' [-dt9 - vA9 + uoM - h • v] , 



6 



FIG. 1: The graphical representation of the propagators of 
the modeL 



where the four terms in the third hne in Eq. (|19p represent 
the Martin-Siggia-Rose action for the stochastic problem 
([3]) at fixed velocity field v, and the first two lines describe 
the Gaussian averaging over v defined by the correlator 
in Eqs. ((10|) and (HI]). 
The action (|19[) is given in a form convenient for a re- 
alization of the field theoretic perturbation analysis with 
the standard Feynman diagrammatic technique. From 
the quadratic part of the action one obtains the matrix 
of bare propagators. The wave-number-frequency repre- 
sentation of, in what follows, important propagators is as 
follows: (a) the bare propagator {99')a defined as 



1 



;o = 



(20) 



and b) the bare propagator for the velocity field {vv)o 
given directly by Eq. pT|) . namely. 



(21) 



where PfjO^) is the transverse projector defined in the 
previous section by Eq. ([5]). Their graphical representa- 
tion is presented in Fig.[TJ 

The triple (interaction) vertex —6'vjdjd = O'vjVjO is 
present in Fig. [21 where the momentum k is flowing into 
the vertex via the auxiliary field 0' . 



IV. UV RENORMALIZATION AND RG 
ANALYSIS 

We start with the analysis of uv divergences, which 
is usually based on the analysis of canonical dimensions 
(see, e.g., [H, [2l,[2l|). First of aU, the dynamical model 
(fT9)) . as well as all models of this type, is a so-called two- 
scale model [IMUlij .29, 30], i.e., the canonical dimension 
of some quantity F is described by two numbers, namely. 




FIG. 2: The interaction vertex of the model (wave-number- 
frequency representation) . 



the momentum dimension d^p and the frequency dimen- 
sion ct^ . To find the dimensions of all quantities it is ap- 
propriate to use the standard normalization conditions 

'k 



d% = di = d'i = 0, 



and the requirement that each term of the action func- 
tional must be separately dimensionless with respect to 
the momentum and frequency dimensions. The total 
canonical dimension dp is then defined ds, dp = dp + 2dp 
[it is related to the fact that dt oc i/gd^ in the free action 
with choice of zero canonical dimension for Uq] . In 
the framework of the theory of renormalization the total 
canonical dimension in dynamical models plays the same 
role as the momentum dimension does in static models. 

The canonical dimensions of the model ([TO)) cannot 
be determined directly because it contains fewer terms 
than fields and parameters. Thus one is faced with some 
kind of uncertainty in calculation of canonical dimen- 
sions. This freedom is demonstrated by the fact that 
parameter ft, = |h| can be eliminated from the action 
(see Ref. [sl] for details). When h is eliminated from the 
action, which is equivalent to the assigning of zero canon- 
ical dimension to it, the canonical dimensions of the other 
quantities can be calculated unambiguously. They are 
present in Table HI where also the canonical dimensions 
of the renormalized parameters are shown. 

The model is logarithmic at e = 77 = (the coupling 
constants go and uq are dimensionless); therefore the uv 
divergences in the correlation functions have the form of 
poles in e, 77, and their linear combinations. 

The quantity that plays a central role in the renormal- 
ization of the model, namely, the role of the formal index 
of the uv divergence, is the total canonical dimension of 
an arbitrary one-particle irreducible correlation (Green's) 
function F = ($ • • • <i>)i_ir- It is given as follows: 



d't 



2d^ 



d + 2- N^dq,, 



(22) 



where iV$ = {N0, Ngr, N^} are the numbers of corre- 
sponding fields entering into the function F, and summa- 
tion over all types of fields is implied. It is well known 
that superficial uv divergences, whose removal requires 
counterterms, can be present only in those Green's func- 
tions F for which the total canonical index dr is a non- 
negative integer. 

A detailed anal ysis of divergences in the problem 
was done in Ref. [S^ (see also Refs. [29, 30] ); therefore 
we shall present here only basic facts and conclusions 
rather than to repeat all details. First of all, every one- 
irreducible Green's function with Ngr < Ng vanishes. On 
the other hand, dimensional analysis based on Table I 
leads to the conclusion that for any d, superficial diver- 
gences can be present only in the one-irreducible Green's 
functions {9'9 ■ ■ ■ 9) with only one field 9' {Ngi = 1) and 
an arbitrary number Ng of fields 9. Therefore, in the 
model under investigation, superficial divergences can be 
found only in the one-particle irreducible function {9' 9). 
To remove them one needs to include into the action func- 
tional a counterterm of the form 9' A9. Its inclusion is 
manifested by the multiplicative renormalization of the 
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TABLE I: Canonical dimensions of the fields and parameters 
of the model under consideration. 



F 


V 


e 


e' 


m, A, 




90 




g,u,h 


d% 


-1 


-1 


d+l 


1 


-2 


2e + v 


V 





dp 


1 











1 











dp 


1 


-1 


d+l 


1 





2e + ri 


V 






bare parameters go,uo, and j/q in the action functional 



= vZ^, go^ g^j^^+'^Zg, uq^uh^'Zu- (23) 



Here the dimensionless parameters g, u, and v are the 
renormahzed counterparts of the corresponding bare 
ones, fj, is the renormahzation mass (a scale setting pa- 
rameter) in the minimal subtraction (MS) scheme, and 
Zi = Zi{g,u) are renormahzation constants. 

The renormahzed action functional has the following 
form: 

S'«($) ^ -\ j dti d'^Xi dt2 d'^X2 

+ J dt d'^x 9' [-dtO - vAe + iyZiA9-h-v], 

where the correlator D^j is written in renormahzed pa- 
rameters (in wave-number-frequency representation) 



pP-{\^)gv^^?''+'^k'^-^-'^^-'^ 



(25) 



By comparison of the renormahzed action (|24p with def- 
initions of the renormahzation constants Z^, i — g,u,v 
[Eq. we come to the relations among them: 

Zv ~ Z\^ Zg = Z^ ^ , Z^ = Z^ ^ . (26) 

The second and third relations are consequences of the 
absence of the renormahzation of the term with D"" in the 
renormahzed action (|24p . Renormahzation of the fields, 
the mass parameter m, and the vector h is not needed, 
i.e., Z^ — 1 for all fields, Z,n ~ 1, and also Zh = 1- 

In what follows, we shall work with two-loop approx- 
imation to be able to see the effects of helicity. The 
calculation of higher-order corrections is more difficult in 
the models with turbulent velocity field with finite corre- 
lation time than in the cases with S correlation in time. 
First of all, one has to calculate more relevant Feynman 
diagrams in the same order of perturbation theory (see 
below). A second and more problematic distinction is 
related to the fact that the diagrams for the finite corre- 
lated case involve two different dispersion laws, namely, 
a; (X fc^ for the scalar field and uj oc k^~'^ for the velocity 
field. This leads to complicated expressions for renormal- 
ization constants even in the simplest (one-loop) approx- 
imation [39, am. But, as was discussed in [33, [39l. 



this difficulty can be avoided by the calculation of all 
renormahzation constants in an arbitrary specific choice 
of the exponents e and rj that guarantees uv finiteness 
of the Feynman diagrams. From the point of calcula- 
tions the most suitable choice is to put 77 = and leave 
e arbitrary. 

Thus, the knowledge of the renormahzation constants 
for the special choice 77 = is sufficient to obtain all im- 
portant quantities like the 7-functions, /3-functions, co- 
ordinates of fixed points, and critical dimensions. 

This possibility is not automatic in general. In the 
model under consideration it is the consequence of an 
analysis which shows that in the MS scheme all the 
needed anomalous dimensions are independent of the ex- 
ponents s and 77 in the two-loop approximation. But 
in the three-loop approximation they can simply appear 

In Ref. [39\ the two-loop corrections to the anomalous 
exponents of model (fT9| without helicity were studied. 
We shall continue those investigations including the ef- 
fects of helicity. 

Now we can continue with renormahzation of the 
model The relation 5(61, 6*', v, eo) = S'^(6l, 6*', v, e, /x), 
where cq stands for the complete set of bare parame- 
ters and e stands for the renormahzed ones, leads to the 
relation M^(A, cq) = W^'{A, e, /i) for the generating func- 
tional of connected Green's functions. By application of 
the operator 2?^ = at fixed cq on both sides of the 
last equation one obtains the basic RG differential equa- 
tion 



(27) 



where I?i?G represents the operation written in the 
renormahzed variables. Its explicit form is 



V 



RG = '^fj. + Pgig,u)dg 



- l3u{g,u)du-"f„{g,u)V^, (28) 



where we denote "Dx = xdx for any variable x and the 
RG functions (the [3 and 7 functions) are given by well- 
known definitions, and in our case, using relations (j26[) 
for renormahzation constants, they have the following 
form 



7i. = V^j^lnZy 



Pg = V^g = g{-2e-ri + Z-f^), 
[in = 'Df,u = u{-ri + j^). 



(29) 
(30) 
(31) 



The renormahzation constant Z^, is determined by the 
requirement that the one-irreducible Green's function 
{6'9)i^ir must be uv finite when is written in renormal- 
ized variables. In our case it means that they have no 
singularities in the limit s,ri — s- 0. The one-irreducible 
Green's function {9'9)i^ir is related to the self-energy 
operator Eg/e by the Dyson equation 



)i-ir = + i^qP^ - ^e'e{uj,p). 



(32) 



Thus Zi, is found from the requirement that the uv di- 
vergences are canceled in Eq. ([5^ after the substitution 
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the second row in Fig.[3]). Here, Sd = 27r'^/2/r(d/2) de- 
notes the d-dimcnsional sphere, 2^1(07 b,c,z) — l + ^z + 

° e'(c+i>i^^'' ^'^ + ■ • • represents the corresponding hyper- 
geometric function. In further investigations the hchcal 
term with in B2 has to be taken with d = 3 but for 
completeness we leave the d dependence in this part of 
B2 in Eq.dSSl). 

Finally, the renormalization constant Zi, = Zi is given 
as follows: 



FIG. 3: The one and two-loop diagrams that contribute to 
the self-energy operator Eg/g. 



i^Q = vZi,. This determines Z^, up to an uv- finite contri- 
bution, which is fixed by the choice of the renormaliza- 
tion scheme. In the MS scheme all the renormalization 
constants have the form (1 -f poles in e, 77 and their lin- 
ear combinations). The self-energy operator Eg'e is rep- 
resented by the corresponding one-irreducible diagrams. 
In contrast to the rapid-change model, where only the 
one-loop diagram exists (it is related to the fact that all 
higher-loop diagrams contain at least one closed loop that 
is built up of only retarded propagators and thus are au- 
tomatically equal to zero) , in the case with finite correla- 
tions in time of the velocity field, higher-order corrections 
are nonzero. In two-loop approximation the self-energy 
operator Sg/g is defined by diagrams that are shown in 
Fig.H 

As was already mentioned, in our calculations we can 
put ?7 = 0. This possibility essentially simplifies the eval- 
uations of all quantities [H, [H, [s^ . Then the singular 
parts of the diagrams in Fig. [3] have the following analyt- 
ical form (for calculational details see Appendix E]) : 



A 



1 / M \ 2£ 1 



(277)^^4^(1-1-^) d 

9 9 



(-) 



(33) 




(34) 



Z^ = l - 



gd-l 1 



1 



e d Au(l + u) £2 ^2 32m2(i + u)3 
f{d-l){d + u) 1 



£ d'^{d + 2) 16u2(l-t-u)5 
X2F1 ( 1,1;2+^; ^ 



„2 3 



P 



2' (1 + u) 

2^1 7:1 TTi TTi 



(36) 



£ 144m2(H-u)3 ^ ' \^2' 2' 2' (1 + u) 



where in the helical part (the last line) we already sub- 
stitute d = 3 and denote g = gSd/i^n)'^. 

Now using the definition of the anomalous dimension 
7i/ in Eq. (P^)) one comes to the following expression: 



7. = -2(5^ + 25^^), 



where 



A 



d-1 1 
d 4u(l-f-u) 



(37) 



(38) 



is the one-loop contribution to the anomalous dimension 
7j/ and the two-loop contribution is 



1 



16d2(d+2)u2(i + „)5 ^ ' ' 2'(l + u)2 

2F^(ih^,JA-x^■ (39) 



144u2(l + u)3 ^2' 2' 2' (l + u)2 



The issues of interest are especially the multiplicatively 
renormalizable equal-time two-point quantities G(r) (see, 
e.g., Ref. [IHl)- Examples of such quantities are the equal- 
time structure functions 



02 



? ? 

g^vp^ 



(27r)2'il6u2(l + u)3 d2 
2F1 (l,l;2 



(d-1) /^ye 1 



d. 1 
2 ' JT+u)^ 



{d + 2){l + u) 



{d - 2)7rp2 



1 1 

2F1 I 2' 2'^ 



2' (l + u)2 



(35) 



Expression A is the result of the one-loop diagram, Bi 
is the result for the first two-loop diagram (the first di- 
agram in the second row in Fig. [3]), and B2 is the result 
for the second two-loop graph (the second diagram in 



Snir) ^ mt,^) ~ 0{t,^')Y 



(40) 



in the inertial range, specified by the inequalities / ~ 
l/A^r^L = l /m {I is an internal length) . Here the 
angular brackets (■ • • ) mean the functional average over 
fields <i> = {0,9',v} with weight expS'^($). The infrared 
(ir) scaling behavior of the function G{r) (for r/l 1 
and any fixed r/ L) 



G(r) ~ zyo^r''<3(r/Z)'^«i?(r/i) 



(41) 



is related to the existence of ir stable fixed points of the 
RG equations (see the next section). In (|¥T|l d^. and c?g 
are corresponding canonical dimensions of the function 
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G, R{r/L) is the so-called scaling function, which cannot 
be determined by the RG equation (see, e.g., Ref. [25j). 
and A(3 is the critical dimension defined as 



lo- 



rn 



Here 7^ is the fixed point value of the anomalous dimen- 
sion 7g = iidf^ln Zg, where Zq is the renormalization 
constant of the multiplicatively renormalizable quantity 
G, i.e., G = ZgG" l36|], and A„ = 2 - 7;: is the critical 
dimension of the frequency with which is defined in 
Eq. (|37p (see also the next section). 

On the other hand, the small-r/L behavior of the scal- 
ing function R{r / L) can be studied using the Wilson 
OPE [25]. It shows that, in the limit r/L — > 0, the func- 
tion R{r/L) can be written in the following asymptotic 
form 



(43) 



10 



-2 
-4. 

-6 
-8 
-10 



FPI 


e=ri/2/ / 

/fpii/ 

FPIV 


FPIII 


(one-loop) 



-10 -8 -6 -4 -2 2 4 



8 10 



where Gf are coefficients regular in r/L. In general, sum- 
mation is implied over certain renormalized composite 
operators F with critical dimensions Ap- In the case 
under consideration the leading contribution is given by 
operators F having the form F„ = {diOdiO)'^ . In Sec. lVIII 
we shall consider them in detail where the complete two- 
loop calculation of the critical dimensions of the compos- 
ite operators F^ will be presented for arbitrary values of 
n, d, u, and p. 



V. FIXED POINTS AND SCALING REGIMES 

Possible scaling regimes of a renormalizable model are 
directly given by the ir stable fixed points of the corre- 
sponding system of RG equations [13, H^- The fixed 
point of the RG equations is defined by /9-functions, 
namely, by the requirement of their vanishing. In our 
model, the coordinates 17* , of the fixed points are found 
from the system of two equations 



(44) 



The beta functions /3g and /3„ are defined in Eqs. ((30)1 and 
(Plj) . To investigate the ir stability of a fixed point it is 
enough to analyze the eigenvalues of the matrix f2 of first 
derivatives: 



d/Sg/dg dPg/du 
d(3u/dg dPu/du 



(45) 



The ir asymptotic behavior is governed by the ir stable 
fixed points, i.e., those for which both eigenvalues are 
positive. 

The possible scaling regimes of the model in one-loop 
approximation were investigated in Ref. [35| . Our first 
question is how the two-loop approximation changes the 
picture of the "phase" diagram of scaling regimes dis- 
cussed in Ref. 35] , and the second one is what restrictions 



FIG. 4: Regions of stability for the fixed points in one-loop 
approximation. The regions of stability for fixed points FPI, 
FPU, and FPIII are exact, i.e., not influenced by loop correc- 
tions. The fixed point FPIV is shown in one-loop approxima- 
tion. The d dependence of the FPIV in two-loop approxima- 
tion is shown in Fig. [5] below. 



on this picture are given by helicity (in the two-loop ap- 
proximation) . The two-loop approximation in the model 
under our consideration without helicity was studied in 
Ref. [3^ but the question of scaling regimes from the two- 
loop approximation point of view was not discussed in 
detail. 

First of all, we shall study the rapid-change limit 
u — > 00. In this regime, it is convenient to make a 
transformation to new variables, namely, w = l/u and 
g' = g/u^, with the corresponding changes in the (3 func- 
tions: 



g'iv 

w{'r/ - 



- 2e- 
7.) 



(46) 
(47) 



In this notation the anomalous dimension 7^/ obtains the 
following form: 

7. = -2(5'^' -t- 2-g'^B'), 

where again g' — g' Sd/ {2tt)'^ . 
A' acquires the form 

A' = -^^^ (49) 



(48) 

The one-loop contribution 



d 4(1 

and the two-loop correction B' is 
(d~l)(dw + l)w2 



B' = 



16d2(d + 2)(H-w)5 



2F, 



w) 



144(1 



■ 2^1 7:1 7:1 77) 



2' 2' 2' (1 



w)2 ^ 

(50) 
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6 

4- 

2 
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-10 
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6=2 

d=3- 

d=4- 



_l 1 I 1 I 1 L. 




-10 -8 



8 10 



FIG. 5: Regions of stability for the fixed point FPIV in the 
two-loop approximation without helicity for different space 
dimensions d. The ir fixed point is stable in the region given 
by the inequalities e > 0, e > t;, and e < d — 1. 



It is evident that in the rapid-change limit w — > (w ^ 
oo) the two-loop contribution B' is equal to zero. This is 
not surprising because in the rapid-change model there 
are no higher-loop corrections to the self-energy operator 
[31L Is^ : thus we are coming to the one-loop result of 
Ref. [35| with the anomalous dimension 7^/ of the form 



7i 



~ lim 



(d-l)g' _ (d-l)g' 



w^a 2d(l -I- w) 



2d 



(51) 



In this regime we have two fixed points denoted as FPI 
and FPU in Ref. [1^. The first fixed point is trivial, 
namely, 



FPI 



0, 



(52) 



with 7^ = 0, and diagonal matrix Q with eigenvalues 
(diagonal elements) 



(53) 



The region of stability is shown in Fig. 21 The second 
point is defined as 



FPU 



0, .9: 



2d 

d-1 



(2e-r;), 



(54) 



with 7* = 2e — rj. These are exact one-loop expressions 
as a result of the nonexistence of the higher-loop cor- 
rections [see discussion below ([50)) ]. That means that 
they have no corrections of order O(e^) and higher [we 
work with the assumption that e ~ 77; therefore it also 



includes corrections of the type 0{rf') and 0{rie)]. The 
corresponding " stability matrix" is triangular with diag- 
onal elements (eigenvalues) 



ni^2{T^-e), = 26 -77. 



(55) 



The region of stability of this fixed point is shown in 
Fig.m 

Now let us analyze the "frozen regime" with frozen 
velocity field, which is mathematically obtained from the 
model under consideration in the limit u — )■ 0. To study 
this transition it is appropriate to change the variable g 
to the new variable g" = g/u [35l|. Then the (3 functions 
are transformed to the following: 



= 5"(-2£ + 27,), 
Pu = u{-ri + j^), 



(56) 
(57) 



where the (3u function is not changed, i.e., it is the same 
as the initial one pip . In this notation the anomalous 
dimension 7^ has the form 



7. - -2i-g"A" + 2-g"^B"), 



(58) 



where, as obvious, g" — g" Sd/ {2'kY . The one-loop part 
A" is now defined as 



J^' = - 



d-l 



1 



d 4(1 + u) 
and the two-loop one B" is given by 
{d~l){d + u) 



(59) 



B" 



2F1 1,1; 2 



16rf2(d + 2)(l + u)5 

-2 /I 1 5 



144(1 + u)3 



2F1 



1 



2 2' 2 (l + u) 



2' (1 + u)^ 
(60) 



In the limit u ^ the functions A" and B" obtain the 
following forms 



A"o 



1 



4d 



(61) 



and 



B"o = 



(d - l)2Fi (1, 1; 2 -f f ; l) np^F^ (i, i; f ; l) 



l6d{d+2) 



144 



(62) 

The system of f3 functions (156]) and ([57]) exhibits two fixed 
points, denoted as FPIII and FPIV in Ref. [s^l, related 
to the corresponding two scaling regimes. One of them 
is trivial, 



FPIII 



(63) 



with 7* = 0. The eigenvalues of the corresponding matrix 
f2, which is diagonal in this case, are 



-2e, ^2 



-V- 



(64) 
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FIG. 6: Regions of the stability for the fixed point FPIV in 
two-loop approximation with helicity. The ir fixed point is 
stable in the region given by the inequalities e > 0, e > r;, and 
e < Ep. 



FIG. 7: Regions of stability for the fixed point FPV in the 
two-loop approximation without helicity. d dependence of the 
stability is shown. 



Thus, this regime is ir stable only if both parameters e 
and r/ are negative simultaneously as can be seen in Fig. ID 
The second, nontrivial, point is 



FPIV 



0, -g': 



e B"o 
'2A"q 2A"l 



(65) 



where A"o and S"o are defined in Eqs. (pT|) and ([S^ . 
respectively. 

First, let us study the influence of the two-loop approx- 
imation on this ir scaling regime without helicity in the 
general d-dimensional case. We denote the corresponding 
fixed point as FFIVq, and its coordinates are 



FPIVo : li* = 0, g'i 



2d 



1 



d~\\ d-l 
with anomalous dimension ^i, defined as 



d-l 

2d 



111 



9* 



9_ 

2d 



(66) 



(67) 



which is the exact one- loop result [35|. The eigenvalues 
of the matrix (taken at the fixed point) are 



ni^2{e 



1 



1-d 



(68) 



The eigenvalue ^2 = = — ry -I- 7* is also an exact 

one-loop result. The conditions g'^ > 0, fii > 0, and 
> for the ir-stable fixed point lead to the following 
restrictions on the values of the parameters e and 77: 



£ > 0, e > 1], e < d — 1. 



(69) 



The region of stability is shown in Fig. [51 The region of 
stability of this ir fixed point increases when the dimen- 
sion of the coordinate space d increases. 

Now turn to the system with helicity. In this case the 
dimension of the space is fixed for d — 3. Thus, our 
starting conditions for a stable ir fixed point of this type 
are obtained from the conditions (|69p with the explicit 
value d=3: e > 0, e > 77, e < 2. But they are valid only if 
helicity is vanishing and could be changed when nonzero 
helicity is present. Let us study this case. When helicity 
is present the fixed point FPIV is given as 



n -II 



16 



(70) 



Therefore, in the helical case, the situation is a little bit 
more complicated as a result of the competition between 
nonhelical and helical terms within two-loop corrections. 
The matrix f2 is triangular with diagonal elements (taken 
already at the fixed point) 



^2 



2e + 



3^ 
16 



(71) 
(72) 



where explicit dependence of the eigenvalue fii on the 
parameter p occures. The requirement to have positive 
values for the parameter g'l, and at the same time for 
the eigenvalues Vli,Vl2 leads to the region of the stable 
fixed point. The results are shown in Fig.[Sl The picture 
is rather complicated due to the very existence of the 
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critical absolute value of p 



(73) 



which is defined from the condition of vanishing of the 
two- loop corrections in Eqs. ([70]) and ([7T|) : 



-1 



SttV 
16 



0. 



(74) 



As was already discussed above, when helicity is not 
present, the system exhibits this type of fixed point (and, 
of course, the corresponding scaling behavior) in the re- 
gion restricted by the inequalities £ > 0, e > 77, and e <2. 
The last condition changes when the helicity is switched 
on. The important feature here is that the two-loop con- 
tributions to g'l and fii have the same structure but op- 
posite sign. This leads to different sources of conditions 
in the cases when \p\ < pc and \p\ > pc, respectively. In 
the situation with \p\ < pc the positiveness of fli plays a 
crucial role and one has the following region of stability 
of the ir-fixed point FPIV: 



e > 0, e > 1], e < 



32 



16 - 37rV2 



(75) 



On the other hand, in the case with \p\ > pc, the prin- 
cipal restriction on the ir-stable regime is yielded by the 
condition > with final ir-stable region defined as 



£ > 0, £ > 77, £ < 



32 



'W + Stt^P^' 



(76) 



Therefore, if we continuously increase the absolute value 
of the helicity parameter p, the region of stability of the 
fixed point defined by the last inequality in Eq. (|75p in- 
creases too. This restriction vanishes completely when 
IpI reaches the critical value pc, and the picture becomes 
the same as in the one- loop approximation [35j. In this 
rather specific situation the two-loop influence on the re- 
gion of stability of the fixed point is exactly zero: the he- 
lical and nonhelical two-loop contributions are canceled 
by each other. Then if the absolute value of the pa- 
rameter p increases further, the last condition appears 
again, namely, the third condition in Eq. (|76p . and the 
restriction becomes stronger when \p\ tends to its maxi- 
mal value, \p\ = 1. In this case of the maximal breaking 
of mirror symmetry (maximal helicity), \p\ — 1, the re- 
gion of the ir stability of the fixed point is defined by the 
inequalities £ > 0,£ = 77, and e < 2.351 (see Fig.lHl). It 
is interesting that the presence of helicity in the system 
leads to the enlargement of the stability region. 

Now let us turn to the most interesting scaling regime 
with a finite value of the fixed point for the variable u. 
But by a short analysis one immediately concludes that 
the system of equations (see also (s^)) 



Pa 



g(-2£ - ry + 37^) 
u{-i] + 7^) = 



0, 



(77) 
(78) 




FIG. 8; Regions of stability for the fixed point FPV in the 
two-loop approximation with helicity in the situation when 
p<Pc = 4/(3^/2^). 



can be satisfied simultaneously for finite values of g,u 
only in the case when the parameter e is equal to 77: 
£ = 7]. In this case, the function f3g is proportional to the 
function As a result we have not one fixed point of 
this type but a curve of fixed points in the g ^ u plane. 
The value of the fixed point for the variable g in the two- 
loop approximation is given as follows (we denote it as 
in Ref. [3| as FPV): 



FPV 



1 

'2A 



1 6* 2 



(79) 



with the exact one-loop result for ^* = e = rj [this is 
already directly given by Eq. ([75| ]. Here and are 
expressions A and B from Eqs. ((38|) and ((39l) which are 
taken at the fixed point value it* of the variable u. The 
possible values of the fixed point for the variable u can 
be restricted (and will be restricted) as we shall discuss 
below. The stability matrix 17 has the following eigenval- 
ues: 



Vii = 0, VL2^ ig* 



dg 



du 



(80) 



The vanishing of f2i is an exact result which is related to 
the degeneracy of the system of Eqs. (ffT]) and (jTS]) when 
nonzero solutions with respect to g and u are assumed, 
or, equivalently, it reflects the existence of a marginal 
direction in the g — u plane along the line of fixed points. 

We start the analysis of the last fixed point with the 
investigation of the influence of the two-loop correction 
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FIG. 9: Regions of the stability for the fixed point FPV in 
two-loop approximation with helicity in the situation when 
P > Pc = 4/(3^/V). 



on the corresponding scaling regime when helicity is not 
present in the system (p — 0). In this situation it is 
interesting to determine the dependence of the scaling 
regime on dimension d. The coordinates of the possible 
fixed points are 



a* 



2du^{d + 



iFi (l, 



1;2- 



2' (l+u,y^ I _2 



(d-l)2(d+2)(l+7/,)2 



£^(81) 



where is arbitrary for now. To have a positive value of 
the fixed point for variables g and u one finds a restriction 
on the parameter e: e > 0. Possible restrictions on the 
ir fixed point value of the variable u can be found from 
the condition > 0. The explicit form of ^i is 



(d-l)(d + 2)(d + 4)(l + u,)6 



X (1 + u*)^(4 + d){2d{u^ - 1) + ("* - 3)-it*) 



X 2F1 1,1; 2 



d 

2' (T 

+4u,(d + u,)2Fi (2, 2; 3 



(82) 



2' (l + u*)2 

In Fig. [71 the regions of stability for the fixed point FPV 
without helicity in the e — u plane for different space 
dimensions d are shown. It is interesting that in the 
two-loop case a nontrivial d dependence of ir stability 
appears, in contrast to the one-loop approximation (35| . 



Now let us turn to the situation with helicity and inves- 
tigate its influence on the stability of the ir fixed point. 
In this case we work in three-dimensional space; thus the 
coordinates of the fixed point are defined by the following 
equation: 



3u*(l + u^)e + 



3M*e 



X (2(3 + u,)2Fi (1,1;^; ^ 



2' {i + u,y 

-5^(l + ^.)V2Fi(i,i;^ 



(83) 



1 



2' 2' 2' (l + u,)2 



The competition between helical and nonhelical terms 
appears again, which will lead to a nontrivial restriction 
for the fixed point values of the variable u to have positive 
fixed values for variable g. Next, the eigenvalue Q.2 of the 
matrix f2 is now 



^2 = 



9 



1 + M* 140(1 
8ii(3 + u)2Fi (2,2-"- 
+14(1 

X 2i^l ( 1, 1 



1 



V 2'(l + «,)2 
^(m*(3 + u*)-6) 
7 1 



2' (1 + w,) 



1 



'3 3 7 



2' 2' 2' (1-f 
1 



2' 2' 2' {l + u,Y 



(84) 



with a nontrivial helical part that plays an important role 
in determination of the region of ir stability of the fixed 
point. 

It cannot be seen immediately from Eqs. ([55]) and ([M]) 
but numerical analysis shows that again an important 
role is played by pc — 4/ (VStt). First let us study the case 
when IpI < pc- The corresponding region of stable ir fixed 
points with g^, > is shown in Fig.[51 In the case when he- 
licity is not present (p = 0, see the corresponding curve in 
Fig. [5]), the only restriction is given by the condition that 
172 > 0; on the other hand, the condition > is sat- 
isfied without restriction on the parameter space. When 
arbitrarily small helicity is present, i.e., p > 0, a restric- 
tion related to positiveness of arises and is stronger 
when \p\ is increasing (the right curve for the concrete 
value of p in Fig. [5]) and comes to play the dominant role. 
At the same time, with increasing of |p| the importance 
of the positiveness of the eigenvalue 0.2 decreases (the 
left curve for the concrete value of p in Fig. [5]). For a 
given IpI < pc there exists an interval of values of the 
variable for which there is no restriction on the value 
of the parameter e. For example, for \p\ — 0.1, it is 
1.128 < < 13.502, for \p\ = 0.5, 0.217 < < 0.394, 
and for \p\ = 0.7, 0.019 < it* < 0.029. Now turn to the 
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case \p\ > pc- When \p\ acquires its critical value pc, the 
ir fixed point is stable for all values of > and e > 0, 
i.e., the condition 0,2 > Q becomes satisfied without re- 
strictions on the parameter space. On the other hand, 
the condition 17* > yields a strong enough restriction 
and it becomes stronger when \p\ tends to its maximal 
value IpI = 1 as it can be seen in Fig.lHl). 

The most important conclusion of our two-loop inves- 
tigation of the model is the fact that the possible re- 
strictions on the regions of stability of ir fixed points 
are "pressed" to the region with rather large values of 
e, namely, e > 2, and do not disturb the regions with 
relatively small e. For example, the Kolniogorov point 
(e = ?7 = 4/3) is not influenced. 

As was already discussed (see the previous section) if 
F denotes some multiplicatively renormalized quantity 
(a parameter, a field, or a composite operator) then its 
critical dimension is given by the expression 



A[F] = A^^ = 4 + ^^d% 



-If: 



(85) 



see, e.g., Refs. |25|,|2 



are the canonical dimensions of F, A, 



3fl| for details. In Eq. ^ d% and 

2 - 7: is 

the critical dimension of frequency, and 7^ is the value 
of the anomalous dimension 7^ = T)^ In Zp at the cor- 
responding fixed point. Because the anomalous dimen- 
sion 7i, is already exact for all fixed points at one-loop 
level, the critical dimensions of frequency uj and of fields 
$ = {v, 6',6''} are also found exactly at one-loop level 
approximation ^35] . In our notation they read 



A^ = 2-2e + ri for FPU, 
A^ = 2-e forFPIV, 
A^ = 2-£ = 2-77 forFPV, 



(86) 



where the values of the critical dimensions correspond 
to the given fixed point (see above in this section and 
TablelU, and / is some scaling function whose explicit 
form is not determined by the RG equation itself. The 
dependence of the scaling functions on the argument r/L 
in the region r/L 1 can be studied using the well- 
known Wilson operator product expansion (also known 
as the short-distance expansion) 24, 25', |2^, [10] • The 
OPE analysis will be studied in Sec. lVIII 



and 



VI. EFFECTIVE DIFFUSIVITY 

One of the interesting objects from the theoretical as 
well as experimental point of view is the so-called effec- 
tive diffusivity 9. In this section let us briefiy investi- 
gate the effective diffusivity P, which replaces the initial 
molecular diffusivity vq in Eq. ([1]) due to the interaction 
of a scalar field 6 with the random velocity field v. The 
molecular diffusivity vq governs exponential damping in 
time of all fluctuations in the system in the lowest ap- 
proximation, which is given by the propagator (response 
function) 

G{t-t',k) = (0(t,k)0'(i',k))o 

= e{t-t')cxp{-i^ok'^{t-t')). (90) 

Analogously, the effective diffusivity D governs exponen- 
tial damping of all fluctuations described by the full re- 
sponse function, which is defined by the Dyson equation 
((32)) . Its explicit expression can be obtained by the RG 
approach. In accordance with general rules of the RG 
(see, e.g., Ref. PH) all principal parameters of the model 
go,UQ, and vq are replaced by their effective (running) 
counterparts, which satisfy the Gell-Mann-Low RG equa- 
tions 



Av = l-7^ Ae = -1, Ae'=d+1. 



(87) 



Now let us consider some equal-time two-point quan- 
tity F(r) that depends on a single distance parameter 
r which is multiplicatively renormalizable {F — ZpF^, 
where Zp is the corresponding renormalization con- 
stant). Then the renormalized function F^ must satisfy 
the RG equation of the form 



{VRG+-fF)F{r)^0, 



with operator I?_rg given explicitly in Eq. (|28p and usu- 
ally "fp = In Zp. The difference between the functions 
F and F^ is only in the normalization, choice of param- 
eters (bare or renormalized), and related to this choice 
the form of the perturbation theory (in Qq or in g) . The 
existence of a nontrivial ir stable fixed point means that 
in the ir asymptotic region r// 3> 1 and any fixed r/L 
the function F{r) takes on the self-similar form 



F(r)^^.^r^-(r/0-^-/(r/L), 



(89) 



dg 



du 



s-r = Pgig,u) , s— = Puig,u 



ds 



ds 



dv _ 

s— = -iy-fu{g,u) , 



ds 



(91) 



(92) 



with initial conditions g\s=i = g,u\s=i = u,i'\s=i = v. 
Here s = k/fj,, (3, and 7 functions are defined in Eqs. pS)) 
- pip and all running parameters clearly depend on 
the variable s. Straightforward integration (at least nu- 
merical) of Eqs. (PT|) gives a method to find their fixed 
points. Instead, one very often solves the set of equa- 
tions /3g((7*,M*) = /3m (5*, It*) = which defines all fixed 
points , It* . This last approach was used above when 
we classified all fixed points. Due to the special form of 
the (3 functions ^U^, ^T^j we are able to solve Eq. ([^^ 
analytically. Using Eqs. (pij) and (|30p one immediately 
rewrites (1921) in the form 



dv 



dg 



2e + 77 — 37,y g 



(93) 
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FIG. 10: The dependence of i/^ on the heUcity parameter p 
for definite ir fixed point values u« of the parameter u. 
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FIG. 11: The dependence of //« on the ir fixed point for 
concrete values of the helicity parameter p. 



which can be easily integrated. Using initial conditions 
the solution acquires the form 



1/3 



gk 



2£ + l) 



1/3 



(94) 



where to obtain the last expression we used the equations 
g^2e+rij^3 _ g^j^3 _ jj^ -y^g emphasize that the above 

solution is exact, i.e., the exponent 2e -I- ry is exact too. 
However, in the infrared region A:<<A~Z~^,g— s- 
(7*, which can be calculated only pertubatively. In the 
two-loop approximation = gi^"^ e + gi^"^ e'^ and after the 



Taylor expansion of g^ in Eq. (|94p we obtain 



\ 1/3 



= 1 - 



(2) 

5* e 



ig 



(1) 



(95) 



Recall that for Kolmogorov values e = rj — A/i the ex- 
ponent in ([95]) becomes equal to —4/3. Let us estimate 
the contribution of helicity to the effective diffusivity in 
the nontrivial point above denoted as FPV ([55)1 . At this 
point e — 7] [(2e -I- ?7)/3 = e] and 



1 



12(1 + u,) 
2(3 + u*) 



^1 ( 1,1; L 



1 



5(l + u*)2" V ' '2' (l + w*)2 
115 1 \\ 



2'2'2' (1-f M*)^ 



(96) 



In Figs.fTOl and fTT] the dependence of the v^, on the he- 
licity parameter p and the ir fixed point u.^, of the param- 
eter u is shown. As one can see from these figures when 

^ oo (the rapid-change model limit) the two-loop cor- 
rections to i^* = 1 vanish. Such behavior is related to the 
fact, which was already stressed in the paper, that within 
the rapid-change model there are no two- and higher-loop 
corrections at all. On the other hand, the largest two- 
loop corrections to v^, are given in the frozen- velocity-field 
limit (m» ^0). It is interesting that for all finite values 
of the parameter there exists a value of the helicity 
parameter p for which the two-loop contributions to 
are canceled. For example, for the frozen-velocity-field 
limit (u* = 0) this situation arises when the helicity pa- 
rameter p is equal to its critical value pc = 4/ (\/37r) (this 
situation can be seen in Fig.fTTj). It is again the result of 
the competition between the nonhelical and helical parts 
of the the two-loop corrections as is shown in Eq. ([96| . A 
further important feature of the expression (|96p is that it 
is linear in the parameter e. Thus, when one varies the 
value of £ the picture is the same as in Figs.[TU| and [11] 
and only the scale of corrections is changed. In Figs.fTOl 
and[TT]we have shown the situation for the most interest- 
ing case when e is equal to its Kolmogorov value, namely, 
e = 4/3. 
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VII. OPERATOR-PRODUCT EXPANSION, 
CRITICAL DIMENSIONS OF COMPOSITE 
OPERATORS, AND ANOMALOUS SCALING 

A. Operator-product expansion 

Let us now study the behavior of the scaHnff func- 
tion in Eq. According to the OPE [H, [H, ill, H^, 
the equal-time product Fi{x')F2{x") of two renormahzed 
composite operators at x = (x' -|-x")/2 = const and 
r = x' — x" ^ can be written in the foUowing form: 

Fi{x')F2{x'') =Y,CF^{r)F.,{^,t), (97) 

i 

where the summation is taken over aU possible renor- 
malized local composite operators Fi allowed by sym- 
metry with definite critical dimensions A^^^ , and the 
functions Cf, arc the corresponding Wilson coefficients 
regular in L^^. The renormahzed correlation function 
{Fi{x')F2{x")) can now be found by averaging Eq. (l97l) 
with the weight expS*^ with S'^ from Eq. ([M]). The 
quantities (Fi) appear on the right-hand side, and their 
asymptotic behavior in the limit — > is then found 
from the corresponding RG equations and has the form 

From the OPE ([W]) one can find that the scaling func- 
tion f{r/L) in the representation ([55)1 for the correlation 
function Fi{x')F2{x") has the form 

f{r/L) = J2AF,ir/L)^-^, (98) 

i 

where the coefficients Ap. are regular in {r/L}'^. 

The principal feature of the turbulence models is the 
existence of operators with negative critical dimensions 
(also known as "dangerous" operators) [H, IH, Is^, |3l|, 
[34| . Their presence in the OPE determines the ir be- 
havior of the scaling functions and leads to their singular 
dependence on L when r/L ^ 0. At this point the tur- 
bulence models are crucially different from the models 
of critical phenomena, where the leading contribution to 
the representation (j98p is given by the simplest opera- 
tor _F = 1 with the dimension Ap = 0, and the other 
operators determine only the corrections that vanish for 
r/L 0. 

If the spectrum of the dimensions A p^ for a given scal- 
ing function is bounded from below, the leading term of 
its behavior for r/L is given by the minimal di- 
mension. As was discussed in Ref. [11], the model under 
consideration belongs to this case for small enough values 
of the exponents e, rj. 

In what follows, we shall concentrate on the equal-time 
structure functions of the scalar field defined as 

5„(r) ^ ([0(x, t) - 0(x')]"), r = |x - x'l, (99) 

which are also interesting from the experimental point of 
view. The representation (j89p is valid with the dimen- 
sions df, = Q and dp = Ap = nAg = —n. In general. 



not only do the operators that are present in the corre- 
sponding Taylor expansion enter into the OPE but also 
all possible operators that admix with them in renormal- 
ization. In the present model the leading contribution 
of the Taylor expansion for the structure functions ([99]) 
is given by the tensor composite operators constructed 
solely of the scalar gradients 

F[n,p] = d,,0---d,J{dm0y, (100) 

where n = p + 21 is the total number of fields 9 entering 
the operator and p is the number of free vector indices. 

B. Composite operators F[n,p]: renormalization 
and critical dimensions 

As the composite operators (|100p play a central role 
in what follows, let us briefly discuss their renormaliza- 
tion. A complete and detailed discussion of the renormal- 
ization of the composite operators is given in Ref. [3^ . 
Therefore, we shall show only the basic moments neces- 
sary to present explicit expressions for composite opera- 
tors. 

The necessity of additional renormalization of the 
composite operators (|100p is related to the fact that 
the coincidence of the field arguments in Green's func- 
tions containing them leads to additional uv divergences. 
These divergences must be removed by special kinds of 
renormalization procedures which can be found, e.g., in 
[2^ [25I . [26| . where their renormalization is studied in 
general. The renormalization of composite o pera tors in 
models of turbulence is discussed in Refs. [28ll30|. Typi- 
cally, the composite operators are mixed under renormal- 
ization. This means that renormahzed operators (which 
are uv finite) are linear combinations of unrenormalized 
ones. In our case, the most important fact in renormal- 
ization of composite operators F[n,p\ is that they mix 
only with each other during the renormalization proce- 
dure; therefore the corresponding matrix of renormal- 
ization constants Z is found from the condition of multi- 
plicative renormalization which can be written as follows: 

F[n,p] = Z[n,p]in',p']FR[n' ,p'] , (101) 

where Ffj denotes the renormahzed counterpart of the 
composite operator F. It is standard to define the matrix 
of corresponding anomalous dimensions as 

7[«,p][n',p'] = ^[7>!p][„",p"]^M^[«",p"][«',p']- (102) 

After the corresponding analysis of diagrams (for de- 
tails see, e.g., Ref. [11]) it can be shown that the renor- 
malization matrix Z[n.p][n' ,p'] i^i Eg. (|10ip is triangular; 
therefore, the matrix of anomalous dimensions (|102p is 
also triangular. Thus, the anomalous dimensions [the 
eigenvalues of the matrix (|102p ] is directly determined 
by the diagonal elements of the matrix (jlOip . namely, 

7[n,p]^'Dt,Z[n,p]in,p]- (103) 
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Our following aim is the calculation of the diagonal 



elements Z\ 



n,p] [n,p] 



of the renormalization constants ma- 



trix Zyn,-p\\n' ,-p'\ ■ H we denote the generating functional of 
the one-irreducible Green's functions with one compos- 
ite operator F\n^j^ [given in Eq. (|100p ] and any number 
of fields d as Y{x\ff) then we are interested in part of 
it, namely, the 0" term of the expansion of r(a;; ff) in 0, 
which will be denoted as r„^p(x;0). Its analytical form 
is the following: 



dxi ■ ■ ■ dxn 0{xi) ■ ■ ■ 6{Xn) 



x{F[n,p]{x)e{xi)---9{xr,)), 
— I dxi- 



dXnO{xi)---9{Xn) 

xr„,p(x;xi, . . . ,a;„), (104) 



where 9{x) is the functional argument, the "classical 
counterpart" of the random field 0. In the zeroth approx- 
imation the functional (|104p coincides with F[n,p] and in 
higher orders the kernel r„^p(x; xi , . . . , x„) is given by the 
sum of diagrams shown in Fig.[T2] (up to two loops). The 
analysis of the diagrams in Fig.[T3shows that for each di- 
agram and for any argument Xi, the corresponding spatial 
derivative can be isolated as an external factor. There- 
fore, using integration by parts, it is appropriate to move 
them onto the corresponding fields 9{xi) in Eq. p04p . As 
a result the functional (|104|) takes the form 



,{x-9) 



1 



dxi - ■ ■ J dxn (xi) - ■ ■ (a;„) 
{x;xi, . . . ,Xn), (105) 



where we define new vector fields ai{x) = di9{x). 

The black circles on the top of the diagrams in Fig.[T2 
represents the vertex of the composite operator F[n,p\. 
Analytically it corresponds to the following expression: 



Vk{x]Xi, ...,Xk) 



5^F[n,p] 



59{xi)...S9{xk) 



(106) 



where k denotes the number of attached lines. It is pos- 
sible to represent it in the following convenient form (see 
Ref. H): 



V{x; xi, . . . ,Xfe) = Jj9ij(5(a; - Xj) 



d''F[n,p] 
Sttji . . . dui, 



(107) 



where is replaced by di9{x) after the differentiation. 

For example, the one-loop diagram shown in Fig. 1121 
has then the following analytical form: 



Ki = J dxi---J dx4:V{x; xi,X2){9{xi)9'{x3))o 

X{9ix2)9'ixi))o{vkix3)viix^))odk9{x3)di9{x4),{108) 

and analogously one can write down all two-loop dia- 
grams in Fig.fT^ 



r(i) 



r(2) = 




FIG. 12: The Feynman diagrams for the function Tn.p in two- 
loop approximation. The Feynman rules are the same as in 
Sec. lIIII The black circle is defined in the text of the present 
section. The one-loop diagram we denote as Ki, the two- 
loop diagrams we denote as K2a (the first diagram) , K2b (the 
second diagram), K2c (the first diagram in the second row), 
K2d (the second diagram in the second row), K2e (the first 
diagram in the third row), and is'2/ (the second diagram in 
the third row). 



To determine the renormalization constants Z, 



it 



n,p] 

is enough to calculate the function T'nlp"^" {x;xi, . . . , Xn) 
with appropriate choice of its arguments Oi since the 
function F^^^ '*" (x; Xi, . . . , Xn) contains only logarithmic 
divergences. A sufficient choice is to replace them by the 
fixed point x, the argument of the operator F[n,p\. Thus 
the expression a^j^ (x) ■ ■ ■ (x) in (jlOSp can be taken out- 
side the integration. As a result we come to the local 
composite operator Tn.p{x;9): 



dxi ' 



(109) 



' dXn ^ [x^ X\ , . . . , Xrfi^. 



After integration one obtains an expression indepen- 
dent of coordinates. The vector indices of F^^p--*'' are 
transformed into combinations of Kronecker 5 symbols 
and their contractions with vector symbols of the prod- 
uct Oil (a;) • • • {x) give the original composite operator 
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F[n,p]. Expression (|108p can then be written, up to an 
uv-finite part, in the form 



Ki = auai 



92 



daidaj 



(110) 



with 

X^^^i = Jdx3 J dMx)0'{x3))odj{d{x)d'{x4))o 
x{vk{x3)vi{xi))o. (Ill) 

Two-loop diagrams fFig. ll2p can be given in the same 
form, namely, 

_ d^F[n,p] ^^^^ 

where z — a,b, c, d. The analytical expressions for X^^' 
{z = a, 6, c, d, e, /) can be easily written in analogy with 
the one-loop diagram Ki ; therefore we shall not give their 
explicit form here. The tensors X^'^'' in Eq. (|112p can be 
decomposed into basic structures made of Kronecker S 
symbols. 



iii2,kik2 



1112,1^1 k2 



K2^ 

111213 jkik^ks 



j ^iii2,kik2> 



(113) 



X 



t4''T^iMk2^ (114) 
2 

X! ^"-^i/i2»3,fcl/C2fc3' (11^) 

3 

j -^iii2iiii,kik2k3kiA-^-^^) 



where z = a,b, c, d, and the tensor structures are defined 
as 



(2) _ SikSji + Sgdjk 

2 ' 



(117) 



T 



(1) 



ijk,lmn ~ g(^^ij{^klSmn + Skm^ln + Skn^lm) 

+ (j^fc) + («^fc)), (118) 



(2) 



ijkjlmn 



(1) 



ijkl,mnop 



-|-(/ to) + (/ <-> n)), (119) 
1 ^ 

+ SmpSon), (120) 



(2) ^ J_ 

ijkl.mnop y2 

(^(^ij (^mn('5fco^;p + (J/tpf^io) + 5mo{5kn&lp + ^kp^ln) 
+ '5„jp(4„(5;o -I- ^o^iri) + <^no('5fcm'5zp -|- Skp^lm) 
+ 5np{5km5lo + Sko^lm) + SopiSkm^ln + Skn^lm)) 

+ [i ^ k) + {j ^ I) + {i ^ k) + {i^ I) 

+ (i^/c,j^/)), (121) 



(3) 



1 



ijkl.mnop 2z|. 
X (^Sim{SjniSkoSlp + SkpSlo) + Sjo{SknSlp + Skp^ln) 

+ 5jp{5kn6lo + Sko^ln)) (122) 

-|- (m ^ n) + (m ^ o) + (to ^ o) + (to <-!■ p)^ , 



and the scalar coefficients A| (x = i^i, K2ai ■ ■ ■ , ^2/) in 
Eqs. (|113p - pi6p are given as 



(d+ l)Xf - 2X^ 
d{d + 2){d-l) 
2{-Xf + dX'i) 
d{d + 2){d-l)' 



(123) 
(124) 



for x = Ki, K2a, K2bT K2C, K2d- Further, for operator 
K2e we have 



K2 



%[d+l)X^^'' ~2X^^') 
d(d + 2)(d-h4)(d- 1) ' 
6(-3Xf^'' + {d + 2)X^^') 
d(d + 2)(d + 4)(d- 1) 



(125) 
(126) 



and for operator K2f 
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A 



A' 



K2f 



9(d + 3) ((d + 5)Xf 



8X 



K2f 



72X 



K2f 



(rf- l)d(d+l)(d + 2)(d + 4)((i + 6) 
72 ({d + ■i)xf'' - + 3d + 6)X^'f + 2{d + 2)xf^'^ 
{d - l)d{d +l){d + 2){d + 4){d + 6) 
24 (3Xf +{d + 2) (-6X^" + id + 4)X^'f'^'^ 
{d - l)d{d +l)(d + 2){d + 4){d + 6) ' 



(127) 
(128) 
(129) 



where Xj , j — l,2,3,x — Ki, K2ai ■ ■ ■ , are given 
in Appendix |D] (calculations are performed in the MS 
scheme). Looking at expressions (|D2I) and (jD16|) in Ap- 
pendix [D] one can see the explicit dependence of the di- 
agrams and K2d on the helicity parameter p. These 
diagrams are not present in the case of the rapid-change 
model at all (m — > cx); see Ref. [s^l) because they con- 
tain closed circuits of retarded propagators {09')o and 
therefore automatically vanish (for the same reason the 
self-energy operator See' of the rapid-change model has 
only one- loop corrections '32'|). This is one of the impor- 
tant reasons to study higher-loop corrections of models 
with finite correlation time, namely, some considerable 
properties cannot be studied within simple rapid-change 
models. 

Let us briefly concentrate our attention on the com- 
parison of the rapid-change limit of our two-loop results 
for composite operators with those obtained in Ref. [s^ . 
This comparison leads to some nontrivial results for the 
corresponding integrals. They are presented in Appendix 
[C] [expressions (|C4)) . ([CS]) . and ([Ce]) ]. We found one mis- 
print in Eq. (5.42) of Ref. [l^l, namely, there must be an 
overall factor c?^ — 1 in the expression for Ai. 

The critical dimensions of our operators are defined by 
th egeneral formula given in Eq. ([55]) . When we rewrite 
it in the concrete form of the operator Fnp = F[n^p\ = 
di^O ■ ■ ■ di^9{di9di9y then we have 



The coefficient C„p will not contribute into the corre- 
sponding anomalous dimension (this can be verified by 
direct calculation); hence we do not present its explicit 
form in what follows. The coefficients Anp and Bnp are 
defined as 



4 — 



(133) 



and 



J^7ip — n-l I 2 1 2 



1,(3) aK2, , 1,(3) 4^2. 



(134) 



where Ay,x = {l, 2},y = {K i,K2z},z = {a, 5, c, d, e} 
are defined in Eqs. (fn3l) - (fT251) and 



= {n~p){d + n + p-2) . 



= n[n — 1) 



(135) 
(136) 



.(3) 

n 

.(3) 



{n-i){n-p){d + n + p~2) , (137) 



k)^' = n{n- l)(n-2). (138) 
Then the anomalous dimensions have the form 



Af„, -4. +A.d5^,. +7f„ . (130) 



Now using the canonical dimensions shown in Table [H 
namely, dg — ~1 and dg=0 one immediately comes to 
the result 



7f„ 



(131) 



This means that the critical dimensions of our operators 
are equal to the corresponding anomalous dimensions at 
a corresponding fixed point. 

The first step to determine the anomalous dimen- 



sions is to calculate the constants Z„ 



[n,p][n,p] 



Eq. (|103p ] in the two-loop approximation. In our case it 



IS given as 



(132) 



Inp = 7F„ 



^2 A. 



np 



9 



(139) 



Thus, the coefficient A^p represents the one-loop contri- 
bution to the anomalous dimension, and the coefficient 
Bnp the two-loop one. The critical dimension l^p^p [see 
Eq. (|13ip ] of the operator Fnp is obtained from (|139p 
when it is taken at the corresponding fixed point. 



C. Anomalous scaling: Two- loop approximation 

Our aim is the investigation of the influence of the he- 
licity on the anomalous scaling in the most interesting 
situation of the degenerate flxed point, namely, the fixed 
point denoted as FPV in SecEl In this case, the di- 
mensions Ai?^ are represented in the following series in 
the only independent exponent £ = [it is obtained from 
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(jl39p by the substitution of the corresponding fixed point 
for g^] 



The one-loop contribution has the form 



(140) 



(1) _ 2n(n- 1) - (n-p)(d + n + p- 2)(rf-H) 
2(rf + 2)(d-l) 

(141) 

which is independent of the parameter u (the ratio of the 
velocity correlation time and the scalar turnover time). 
Although the fixed point value given by Eq. (j83p and 



the coefficient Bnp in Eq. (|139p explicitly depend on the 



helicity parameter p, the two-loop contribution to the 
critical dimension A]^ is independent of p. Thus, the 
result is the same as that obtained in Ref. [s^ (there is 
a misprint in the final explicit result but the correct for- 
mula was republished in Ref. [76|l. Its explicit expression 
is rather large, and as it can be found elsewhere we shall 
not repeat it here. At first sight this result is a sur- 
prise but it can be understood in the following, rather 
simple, way. As we know the structure functions S'n(r) 
(which are studied here) are functions of the value of the 
distance r = |x — x'|. Therefore, only those phenom- 
ena will have impact on the critical dimensions that can 
"change" the spatial distances. Among such phenomena 
belong the compressibility and anisotropy. As for helic- 
ity, it breaks the mirror symmetry but it does not disturb 
spatial distances. Therefore, it cannot influence the crit- 
ical dimensions, i.e., it cannot change the corresponding 
asymptotic behavior. Thus, if our statement is right then 
we expect that the situation will be the same in all orders 
of perturbation expansion, namely, the quantities such 
as effective diffusivity will depend on helicity, but criti- 
cal dimensions of the structure functions will not. But, 
of course, for now it is only a speculation and the inde- 
pendence of the critical dimensions of helicity is maybe 
only the effect of two-loop approximation. To solve this 
problem at least three-loop calculations are needed. 

On the other hand, to study the helicity effect on the 
two-loop level it is enough to avoid the conditions of 
isotropy or incompressibility of the system. Thus, the 
next step is, e.g., to include the assumption of compress- 
ibility of the system and investigate the combined effects 
of the helicity and compressibility on the scaling proper- 
ties of the model under consideration. We assume that a 
nontrivial result can be obtained. 

As was already mentioned detailed analysis of the two- 
loop contribution to the critical dimensions of the struc- 
ture functions within the model under our consideration 
(without helicity) was done in Ref. ^SOj]. We have recal- 
culated their results and found some discrepancy in in- 
terpretations of our and their numerical results. That is, 
as our calculations show, a hierarchical behavior of the 



quantity 



/•n? as a function of n 



for a concrete value of d (dimension of space) is destroyed 
in Fig. 1 (d) in Ref. . This figure corresponds to a large 




FIG. 13: Behavior of the quantity C,n (see text) for n = 
4, 6, 8, 10, 20, and 50 as a function of u for d = 10 in units of 



10" 



In Ref. [39| the curves are interpreted as n = 4, 6, 8, 20 



(from bottom to top). But in fact they correspond to n = 
20, 4, 8, 6 (from bottom to top). The curves for n = 10, 50 are 
added to demonstrate the situation more completely. 



enough value of d (namely, d — 10). Our calculations lead 
to the same curves as theirs but they correspond to dif- 
ferent values of n which can be seen by direct comparison 
of Fig. 1(d) in Ref. [3§| and Fig. [13] in the present paper. 

The conclusion is as follows: The hierarchical behavior 
of the quantities Cn is not present for large enough space 
dimensions d (the same situation also occurs for other 
large values of d as can be shown by direct calculations). 

All the other results are the same as in Ref. [39J ; thus 
we shall not repeat them here. 



VIII. CONCLUSION 

In this paper, the infiuence of helicity on the stability 
of asymptotic regimes, on the anomalous scaling, and on 
the effective diffusivity was studied in the framework of 
the passive scalar advected by the turbulent flow with 
finite correlations in time of the velocity fleld. Such in- 
vestigation is important and useful for understanding of 
efficiency of toy models (like the Kraichnan and related 
model) to study the real turbulent motions by means 
of modern theoretical methods including renormalization 
group approach. Thus, it can be considered as the first 
step in investigation of the influence of helicity in a real 
turbulent environment. 

In the present paper, the RG calculations are per- 
formed in the two-loop approximation, which is neces- 
sary to include the effects of helicity. It was shown that 
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the anomalous scahng of the structure functions, which 
is typical for the Kraichnan model and its numerous ex- 
tensions (see the Introduction), is not changed by the 
inclusion of helicity in the incompressible fluid. This is 
given mathematically by the very interesting fact that 
although separated two-loop Feynman diagrams of the 
corresponding composite operators strongly depend on 
the helicity parameter p, their sum - the critical dimen- 
sion A„ - is independent of p in the asymptotic regime 
defined by an ir stable fixed point. This very interesting 
fact can be explained physically by rather simple argu- 
ments in the following way (as was discussed in the pre- 
vious section maybe it is only the effect of the two-loop 
approximation, therefore, to confirm what follows higher- 
loop calculations are needed). The single-time structure 
functions Sn (f ) of the scalar field depend only on the spa- 
tial distance r = |x — x'| but not on the direction. Thus, 
we suppose that only this phenomenon will change the 
critical dimensions of structure functions which modify 
spatial relations. It can be, e.g., the inclusion of com- 
pressibility or spatial anisotropy. On the other hand, 
helicity breaks the mirror symmetry, which is not re- 
lated to distance. As a result, the critical dimensions 
of the structure functions are not affected by the helicity 
of the system. We suppose that an analogous situation 
will hold for all quantities that depend only on the spa- 
tial distance (more precisely, that are constructed from 
quantities taken at different spatial points), and in all 
orders of perturbation theory. The situation can be dif- 
ferent when one includes in the investigation the effects 
of helicity together with another assumption about the 
turbulent flow, e.g., its compressibility. We suppose that 
nontrivial results can be obtained in this more general 
case. 

On the other hand, the stability of possible asymp- 
totic regimes, the values of the fixed RG points, and the 
turbulent diffusivity strongly depend on the amount of 
helicity. It is shown that the presence of helicity in the 
system leads to restrictions on the possible values of the 
parameters of the model. The most interesting fact is the 
existence of a critical value pc of the helicity parameter p 
which divides the interval of possible absolute values of p 
into two parts with completely different behavior. It is re- 
lated to the existence of a competition between nonhelical 
and helical contributions within the two-loop approxima- 
tion. As a result of this competition, within the so-called 
frozen limit, the presence of helicity enlarges the region 
of parameter space with a stable scaling regime, and if 
\p\ = Pc the corresponding two- loop restriction vanishes 
completely and one comes to the one-loop results [35| . 
Similar splitting, although more complicated, into two 
nontrivial behaviors of the fixed point was also obtained 
in the general case with finite correlations in time of the 
velocity field. 

Another quantity which rather strongly depends on the 
helicity parameter p is the effective diffusivity. It is shown 
that the value of effective diffusivity can be 50% larger in 
the helical case in comparison with the nonhelical case. 
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APPENDIX A 

In principle a few ways exist to evaluate the two-loop 
diagrams Bi and B2 which are explicitly shown in Fig. [31 
We compare two of them on the example of the Feynman 
diagram Bi . The explicit analytical expression for Bi in 
the wave-number-frequency representation is 



d'^kd'^qdLLikdLLiqk 



i-d-2eq1-d-2e 



^1 (27r)2rf+2 j (tj2+j,2y2/.4)(^2_^j,2y2^4) 

P»l {P ~ fc)»3 jp-k- g).,4 jP - k)i^P[^i^ (fc)-^3»4 (g) 

■ {iujk + v{v - \^YY{i{uJk +ujq) H- i^(p - k - q)2) 



(Al) 



where p denotes the external wave number (momentum) , 
k and q are two independent internal wave numbers, the 
external frequency ojp is taken to be zero (we are only 
interested in the divergent part of the diagram and it 
is independent of LOp), uJk and ujq are the corresponding 
internal frequencies, 77 is taken to be zero (see discussion 
in Sec. lIVp . P[j is the helical transverse projector defined 
in Secini and over the internal vector indices — 
1, 2, 3, 4, corresponding summations are assumed. 

After integration over the internal frequencies, which is 
rather simple, and then making an expansion in respect 
of the external momentum p and leaving only terms of 
order two in respect of p (because uv divergences can 
have only this structure; see Sec. lIV|) . and, in the end, 
after summation over vector indices one has 



Bi 



-2e 



2 2 
P 1 



(27r)2'i4u2(H-u)2 J ((l + u)(fc2 + 52)4.2k.q) 
(p-k)2q2 p2(]^.q)2 (p.k)2(k-q)2" 



.(A2) 



Now we shall calculate Bi in two different ways. 

The first approach to the integral Bi is based on di- 
vision of integrations into radial and angle parts which 
after standard procedures leads to (for details see, e.g., 
Ref. [13) 

^ g'^vp'^^p^SdSd^i d ~ 1 
^ (27r)2rf4u2(l + u)2 d ^ ' 



dk 

fcl+25 



dq 

3+2e 



L (i„^2)(d-l)/2^2 

fiT 

1 (H-M)(P+g2)+2fcgx' 

where x = cos '0, and tp denotes the angle between vectors 
k and q. Sd and 5^-1 are d-dimensional and {d — 1)- 
dimensional sphere, respectively, which are defined as 
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Sd = 2TT'^/^/r{d/2). The needed ir cutoff of the integra- 
tions is represented by m ~ L~^. It is useful to rewrite 
the denominator of the expression under the integration 
in the form of the infinite series 



(1 + u){k^ + g2) + 2kqx ^ [(P + q^){l + u)]i 



{-2kqx) 



(A4) 



Now we are able to integrate over the angle variable x 
cos term by term in the sum, which yields 



g^^MV^d^d-i d-l fd+l 



dk r°° dq q^{2kqy-^ 



(A5) 



The factor [(—1)^^^ + 1] keeps only odd terms of the 
series. Therefore, we can redefine the summation in 
Eq. (|A5|) in the following form: 



92(2A;g) 



2j 



-■ (A6) 



dk 

To proceed it is appropriate to make transformation of 
the variables k, q to the new polar coordinates, namely: 
k = s cos (j),q ~ s sin (j). This gives 



(27r)2'i4u2(l + u)2d 

g r(2i±i) 22^- 



rT/2 



p l^d+22±2'^ (1 + u)2j + l J^^^ sl+4e 

# cos 02^-1-2^ sin 02J+i-2e, (A7) 







The integral over the radial variable s is trivial, and the 
integral over angle (j) can be found, e.g., in Ref. [73|. Then 
one has 



Bi = 



g r(2i±i) 



(27r)2'i32u2(l + M)2d 

22- rO'-e) 



\niJ e 



(A8) 



In the end, the summation over j leads to the final result 
for Bi which is given in Eq. , where only the divergent 
part is shown. F (j — e) in Eq. (jA8|) for j = has a pole 
with respect to e, which gives the pole of the second order 
in Bi in Eq. 

The second approach to the calculation of our two-loop 
diagrams is as follows. We start with the expression for 



Bi as represented in Eq. I|A2|) . Now using the well-known 
formula of Feynman parametrization 



4"i 4"2 ... Aa„ 



dui ■ ■ ■ dun- 



n 

nr(a,) 



i=l 

5 



1) n< 

i=l 



^ 1=1 



(A9) 



Bi obtains the form 



B, 



(27r)2'i4M2(l + u)2 7q 



duiul 



rfkdq ( p2 — 



(p-k)' 



(k-q)^ 



where 



/c''+2-[Xg2 + 2Yk ■ q + Zfc2] 2+'^+i ' 



X = ui + (l-ui)(l-t-u), 
Y = 1-ui, 
Z = + 



(AlO) 



(All) 



The integration over q is now done using the general 
formula (|B1[) given in AppendixlBl namely. 



dq 



[Xq2 + 2yk-q+Zfc2]l+'^+i F(f 

1 r % F(£ + i)y2 fc,,fc/ 
2 {xz - y2) fc2 



k^^z-^y 

which yields 



+ 1) 
(A12) 



Bi = 



g^v^i^^TT^jd^l) T{e) 
(2^)2'i8w2(l+u)2r(|+e) 



dui—^ 



X-2- 



[z-^-i) 



(A13) 



The integration over vector k is now straightforward, af- 
ter which we have 



Bi = 



A*\4eF(e) 



(27r)2'i32u2(i + u)2rfr(f +£) Vm/ e 

£-1 



du\ — 



X-2- 



(AM) 



where m ^ represents the needed ir regularization. 
We are interested only in the divergent (poles in e) part of 
the expression (jA14p . After doing all needed expansions 
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with respect to s one has the final resuh for the Feynman 
diagram Bi: 



Bi 
1 

f- 

e 



(27r)2d32u2(i + u)3rf2 
^(1) - (d/2) - 



d{l + u) 



■-1 



(m) { 



2' 1 



1 w{ 
dui — 



In Z 



X 



(A15) 



where ^{n) is the logarithmic derivative of the gamma 
function, given by 'ip{n) — T' (n) /T{n). 

Thus, we have two different analytical representations 
of the same Feynman diagram namely, Eqs. ([55)1 and 
(jA15|) . The comparison of them leads to the nontrivial 
result for integral in Eq. (jAlSp . which is present in Ap- 



pendix [C] in Eq. (jCip . The comparison of the results for 
diagram B2 (see Fig.[31) obtained by these two methods 
yields other nontrivial results for some integrals. They 
are shown in Appendix [Cl in Eqs. (|C2p and (jC3[) . 



APPENDIX B 

In this appendix we introduce the general integral for- 
mula which was used in the previous appendix. 

Theorem: Let V be a d-dimensional Euclidean vector 
space over the field of real numbers M. Let ?,n G N 
(natural numbers), and k*^'^ for i = 1, 2, • • • , Z are vectors 
in V. Then for an arbitrary I x I real matrix Vjs with 
detu 7^ 0, arbitrary vectors a^'^ (i = 1,2,- •• ,?), and 
arbitrary c, a € M the following general formula holds: 



[wi.kW -kl^) -|-2aW -kW + c]" ~ 
(-l)"7rf (deti.)-* [}i r(a-f -LtJ+p)[c-(i;-i),saW-a('^)]'^'^^ 



00 J —00 



LfJ + f-"-p 



Via) 



E 

p=0 



([fj - py. i2p + n{mod 2))!4LfJ-P 



^ («-^),...a(-)(.-i),,.,ai-'...(.-i) 
P(qiii ■.q2i2,---,qnin) 



('**2p + n(mod 2)) 



92p + n(mod 2)S2p + Ti(mod 2)^2 



2p + n(mod 2) 



^ '^i2p + 7i(mod 2) + li2p + n(mod 2) + 2 ) 



'j2p + n(mod 2) + 1 92p + n (mod 2) + 2 " ' 



(BI) 



, in , k'f^ and a^*'' are the jth compo- 
5ij denotes the Kro- 



where the summation is taken over all permutations of 
the indices ^l,^2^ ■ 
nents of the vectors k^^^ and a^* 
necker delta, and [n/2j — n/2 for even n and [n/2j 
(n — l)/2 for odd n. Over all dummy indices the corre- 
sponding summation is assumed. 

We shall not present a detailed proof here because it is 
rather large although straightforward; instead we give a 
short recipe for it. To prove formula (|Bip it is appropri- 
ate to use mathematical induction. First, the theorem 
is correct in the scalar case [n = 0). In this specific 
situation the formula is well known (see, e.g., Ref. (25i] ) 



bersome algebraic manipulations the formula for n 
obtained. 



1 is 
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dk(i) . . . dk(') 



[i;,^kW •k(^) +2aW -kW +c]" 



Trf (detz;)-ir(a - f) p 



{v- 



,(■0 . 



(B2) 



Now, let us suppose that formula (|B1[) is valid for n e 
N, n > 1. Then if one differentiates both sides of 



Eq. (jBl[) with respect to a 



(g^+i) 



together with some cum- 



in this appendix we present the integrals that were 
obtained during the calculations of two-loop Feynman 
diagrams Bi and B2 which are shown in Fig. [3] (see Ap- 
pendix]^ for details of the calculations). They are given 
in Eqs. (|Cip - (|C3|) . In addition we present here the an- 
alytical expressions for integrals that were obtained by 
comparison of our two-loop results for composite opera- 
tors with those obtained within the rapid-change model 
[13 [Eas. (fC4ll - (IC6|l ]. 
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1 hn({l + u){l-ui) 
dui 



«i+(l+u)(l-ui) 







dil + u)\ 



f y{l)-^P{d/2) 

2Fi(l,l;2+|;^yqi^ 
(d + 2)(l + w)2 

[u(2 + m) - + u - 1)] 



(Cl) 



dui 



id + 2){i + uy 



(C2) 



1 , Ui[{ui + (1 + W)(l - Ui))(l + U) - (1 - Ml)] 



(Ui + (1+U)(l-Mi))2(l- 7.1)5 



1 1 5 



3(l + u)' 'V2' 2' 2' (1 + M)2y' 



(C3) 



dx (1 - a;2)'*/2 a; arctan 



f l + x 



2 Sd-i d{d + 2)2 ' 



n _^2)(d+i)/2^ 

aa; ^=^= arctan 



/-I V4 - .x2 

1 5d 



2 + a; 



SS'd-i d(d + 2)(d + 4)- 



i^i(l,l;3 



a; 
d 1 
2' 4 



(C5) 



1 (i_^2)(<i+l)/2^3 



da; 



1 



\/4 - 



arctan 



/ 2 + x 



V\/4 - x2 



8 5d-id(d + 2)(d + 4) 



Fi 1,1;3+-;- - 



2'4/ d + 6 



^^i(i,l;4+^;i 
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The explicit form of the coefficients X'^ , j = 1, 2, 3, a; = 
Ki,K2a, ■ ■ ■ ,K2f (see Sec lVII|) . in the MS scheme is the 
(C4) following: 



Sd d-1 g ^K,^Q 

lT:Y\m) 4M(l + u)e' ^ 



.9 '^d 



vK2a ^ 

1 (27r)2'il6u2(l + M)(l-u2) d 



(27r) 

d - 11 ( 



(Dl) 



i^i i,i;2 



d. 1 



2 ' (1+n)^ 



M2^^i (i,i;: 



d. 1 



2' 2(l+it) 



(d + 2)(l + w) 



1 



2'2' 2'(l + ?/)2; ^2{\ + u) V2' 2' 2'2(l + u) 



11^ d 



-R'2a 

2a 



n2<?2 
9 



V^2a V^2a i V^2a 

^2 — ^2a + ^2b ' 



(d- l)(d+ 1) 1 ^/i^^^ 

771/ 



(27r)2rf327i3(l+?l)(l-7l2) d(d+2) £ 



3i^2 7^,1,1; 7^,2 



1 



X 



g^SdSd-1 



-(-) 



4e 



l\-3F2[ ^,1,1; ^,2 



2' (1 + 77)^ 



(27r)2'il67i(l + 7i)(l - 7*2) e V777/ 7_i (1 - 7i)2 + 47tx2 

^ 2(1 + 7.)x arctan (^) 2(1 + 377)a: arctan (^-^=^^=5) 



v/2(l + 7i) - a;2 



+ (77 - 1) In 



1 + 71 



(D2) 
(D3) 

(D4) 



(D5) 



25 



■^1 — ^la + ^Ib 



(D6) 



X 



22Fi(l,l;2+|;^) „ln(i±H)~ 



(27r)2'i32u2(l +w)2 d \mJ e\e {d + 2){1 + u^il - u) (1 - w) 
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